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ABSTRACT 

Based on the Solar Standard Model SSM of Bahcall and Pinsonneault (SSM-BP2000) we developed a solar 
model in hydrostatic equilibrium using three polytropes, each one associated to the nuclear, the radiative and 
the convective regions of the solar interior. Then, we apply small periodic and adiabatic perturbations on this 
tri-polytropic model in order to obtain proper frequencies and proper functions which are in the p-modes range 
of low order < / < 5; for I = 2, 3 and 4 these values agrees with GOLF observational data within a few percent. 

Subject headings: Solar Standard Model (SSM), Lane-Emden, Non Radial Oscillations (NRO), p-modes 



1. Introduction 

Polytropic models have largely been used in the 
study of NRO of a gaseous sphere [Cowling, 1941; 
Kopal,1949; Scunaire,1974; Tassoul,1980]. 

We have computed the first modes of a tri- 
polytropic model (TPM) whose indices n\ = 1.50, 
ri2 = 3.78 and 113 = 20 describes the convective, 
radiative and nuclear zones respectively of the 
solar interior. We used Cowling's approximation 
[Cowling, 1941] which reduces the order of the sys- 
tem of differential equations to 2 (instead of 4). 
The radial part of the perturbation obeys equa- 
tions (1) and (2) [Ledoux-Walraven,1958]: 
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are proper functions, I is the degree of the spherical 
harmonic, Y\ is the adiabatic exponentent equal 
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to I, a is the angular frequency, N is the Brunt- 
Vaisala frequency and P; is the Lamb frequency. 
The equations (1) and (2) and the boundary condi- 
tions lead to a eigenvalue problem with eigenvalue 
cr 2 , which is the problem to be solved. 

2. Polytropes 

Our unperturbed model consist of a gas of par- 
ticles with spherical symmetry, selfgravitating, in 
hydrostatic equilibrium and with its state equa- 
tion given by: 



P = Kp 1 = Kp 
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K and 7 are parameters that depend only on the 
polytropic index n, and the mass and radius of the 
configuration. The polytrope theory developed by 
the ends of the XIX century, can be used to know 
the dynamical structure of a star, within which 
local quasistatic thermodynamic changes follows 
a polytropic process, i.e. one in which the spe- 
cific heat remains constant. This approach can be 
used in some regions of Sun's interior. We have 
used the pressure and the density data from the 
sophisticated SSM of Bahcall and Pinsonneault 
to plot 7 = d Jf lP vs x with x = t/Rq. Vari- 



ous regions clearly emerge. Of these regions, the 
outermost one (71 = 1.677, i.e. n\ — 1.50) repre- 
sents the convective zone where heat transport is 
achieved by adiabatic convection. It's SSM output 
in Fig. 1 is approximated rather well by a constant 



straight line indicating a polytropic behavior. The 
second intermediate zone labeled radiative zone in 
Fig.l can be approach by a polytrope 72 = 1.264, 
i.e. ri2 = 3.788. Finally, we have taken an ave- 
rage of 7 for the innermost regions labeled nuclear 
zone, so we represent this zone by a polytrope with 
73 = 1.050, i.c «3 = 20.0. This three polytropcs 
have been use by Hendry (1993). 
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Fig. 2. — Solutions of the Lane-Emden equation. The 
center of the Sun as at V = 0, U = 3 and its surface at 
U = 0, V --= 00. We choose 9' c _ = -0.00325 (curve C) as 



. We choose 9' e 
a M-Solution for the convective zone and <j>' = 
(curve E) as a M-Solution for the radiative zone. 
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Fig. 1. — 7 at the solar interior due to SSM [Bahcall and 
Pinsonncault,2000]. Regions with 7 constant are described 
by polytropic process. The tri-polytropic model is obtained 
if one assume three values for 7; | in the convective zone, 
1.264 in the radiative zone and 1.050 in the nuclear zone. 

2.1. Three polytropes within the Sun 

Following Hendry [Hendry, 1993] we use £, 9 as 
the variables in the Lane-Emden equation for the 
convective zone with index m; i], 4> as the variables 
for the radiative zone with index n 2 and £, ip as 
the variables for the nuclear zone with index n%. 
Hence the parametric polytropes are: 
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The main challenge is to learn how to fit these 
three polytropes together. Since the physical 
quantities P, p and M are continuous across the 
interfaces (not for example 9 or <j>), the variables 
U and V given by 



V = 



(n + in-e 1 ) 
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Let us start by considering the convective zone. 
We take m = 1.5 thus £1 = 3.6538 and 0^ = 
—0.2033. Though, this polytrope is not being used 
in the vecinity of £ = 0, it is possible to consider all 
of the solutions of the Lane-Emden equation for 
n\ = 1.50. These may be generated beginning at 
£1 with an arbitrary starting slope and integrating 
inwards. Solutions with starting slopes less nega- 
tive than 91 are of particular interest (in the lite- 
rature, they are referred to as M-solutions) since 
these are the ones which intersect the polytrope 
that represents the radiative zone n 2 = 3.788. 
Four such solutions, translated into the U, V vari- 
ables, are shown in Fig. 2, as the curves A, B, C 
and D there. Solutions with starting slopes more 
negative than 0£ (F-solutions) do not intersect 
the radiative polytrope and so do not need to be 
considered here. In the same way we have worked 
the intersection between the nuclear and radia- 
tive zones; we have fixed the nuclear polytrope 
with the #20 E-solution which intersects the M- 
solutions (Fig. 2, as curves E and F) that repre- 
sents the radiative zone n 2 — 3.788, this occur at 
the neighbors to 0.25i?Q. 

Knowing 6*(£), <j){rj) and VKC) we can deduce the 
density and pressure curves p(r) and P(r). The 
above model yields a central density 3 of p c = 
1.483 x 10 5 Kgm in comparation with p c = 
1.299 x 10 5 Kgm~ 3 obtained for the bi-polytropic 



are be very useful [Chandrasekhar,1939]. 



3 SSM value: 
neault,2000] 



1.524 X 10 5 Kgm 3 [Bahcall and Pinson- 



model BPM [Pinz6n-Calvo-Mozo,2001]. 
A comparison of our BPM and TPM models 
with the SSM-BP2000 arc shown in Fig. 3, show- 
ing a better fit of the TPM with respect to the 
SSM-BP2000 in the most central part. 
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Similarly, the value of the Lamb frequency associ- 
ated to the mode I is given by: 
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Fig. 3. — Density (in g X cm 3 ) at the solar interior 
for TPM (thin line), for the BPM (dashed line), for the 
SSM-BP2000 (thick line) and for the n = 3 one polytrope 
model. 



2.2. Characteristic Frequencies 

The characteristic N and L frequencies so called 
Brunt- Vdisald and Lamb frequencies associated 
with restoring forces of pressure and gravity in a 
polytrope of order n can be calculated as functions 
of radius in terms of the Lane-Emdcn function 6. 
Thus, if we define: 






(11) 



where Rq and Mq are the sun mass and the sun 
radius, and G is the gravitational constant, we can 
write a normalized Brunt-Vdisdld frequency as: 



L 2 
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where £ is the Lane-Emden coordinate. 

3. NRO in a Tri-polytropic Model (TPM) 

Space oscillation properties of the solutions of 
equations (1) and (2) are related to the signs of the 
coefficients given by the second members of these 
equations. Space oscillations are allowed only in 
the regions where these coefficients have opposite 
signs. The limits of these regions arc defined by 
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where CC is the central condensation of the poly- 
trope, defined as the ratio of central density to 
mean density. The parameter no is the effec- 
tive polytropic index associated to the pulsations 
[Mullan-Uhich,1988], which in turns related to the 
adiabatic exponent 4 T by: 
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Fig. 4. — Propagation Diagram (I = 2) for the BPM 
(thin lines), for the polytropic index 3 model (thick lines) 
and for the SSM-BP2000 (dashed lines). The TPM 
agree exactly to the BPM except near the center; in this 
region some points for the TPM are shown. 



In the (x,u> 2 ) plane, these equations define two 
curves. In Fig. 4 we have plotted (15) and (16) as 
thin black curves for the BPM. We can see that 
the Lamb frequency diverges near to the center; 
conversely, the Brunt- Vdisald runs from zero at 



the origin and diverges near the surface. This fea- 
ture is common in all models. The TPM propa- 
gation diagram is very similar to BPM except 
at the nuclear region, where there is a prominent 
change. We have been denoted p-modes and g- 
modes the regions of (x,cl> 2 ) plane corresponding 
to the conditions of position and frequency in the 
star, allowing spatial oscillations. These regions 
are characterized by the possibility of existence 
of progressive acoustic waves and progressive gra- 
vity waves respectively [Scuflaire,1974]. Thus we 
shall refer to these regions as the acoustic and the 
gravity regions. We have also plotted in the same 
figure the frequencies (horizontal dot lines) for the 
first p-modes from pl2 — mode for the TPM. In 
addition in the same figure, we can see the propa- 
gation diagram for the SSM [Bahcall and Pinson- 
neault,2000]. 
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The regularity condition at the centre, first re- 
quires that 



u> y 



lz = 



(24) 



and second, the cancellation at the surface of the 
lagrangian perturbation of the pressure which can 
be written down as: 



-y- Z = 



(25) 




Fig. 5. — Behavior of the boundary condition at surface 
(equation (25)) for TPM model as a function of a in /iHz. 



The equations (1) and (2) are very convenient for 
the analytical discussion, but for the numerical 
computations we use the more appropriate form: 
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- = -^y-l z]+Re A^y-z) (18) 
where we have put: 
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In order to determine the solution uniquely, we 
impose the normalizing condition 



y 



l 



(26) 



at the centre. With a trial value for uj 1 we 
integrate equations (17) and (18), with initial 
conditions (24) and (26) using Runge-Kutta 
method, with a step size taken from a pa- 
per of Christcnscn-Dalsgaard, equation (A55) 
[Christensen-Dalsgaard and Mullan,1994]. Usu- 
ally this solution does not satisfy equation (25) 
and a new integration is performed with another 
value of u> 2 . This procedure is repeated until equa- 
tion (25) is satisfied, using a Newton-Raphson 
method to improve the value of lo 2 . 



3.1. Phase Diagram 

The radial displacement S(r) and the pressure per- 
turbation P' are periodic space functions; the vari- 
ables v(r) and w(r) vary strongly from the center 
to the surface, then it is impossible to plot them 
directly along the axes. However, the most appro- 
piate functions: 
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Fig. 6. — Some p-modes for polytropic models. 

Figures A and B corresponds to pi2 modes with a„ y i = 
012,2 = 1796.89 /iHz and 012,4 = 1880.45 /iHz respectively 
for the n = 3 one polytrope model. C (012,2 = 2059.86 
fJ.Hz), and D(oi2,4 = 2165.79 fJ.Hz), are the same modes 
but in BPM, while E (021,2 = 3368.89 fiHz) and F 
(021,5 = 3531.43 fiHz) represents the P21 mode for TPM. 

have been plotted (27) and (28) in each axis of 
the Fig. 6; their signs are chosen according to the 
signs of the variables 6(r) and P'. The number 
of intersections of each curve with the ordinate 
axis in Fig. 6 (the origin excluded) is equal to the 
order of the mode. The sense of rotation of this 
curves is a characteristic feature for the p-modes 
[Scuflaire,1974]. 

4. Conclusion 

Although we do not use an atmosphere model 
and the input physics is described by a multi- 
polytropic structure, the modes obtained from 
the TPM are close to the observacional data. 
We show in Fig. 7 a comparison with GOLF 
data 5 , showing an agreement within a few percent 
(roughly between 4.5% and 7%) up to a radial 
order of 30. However, we can see that the TPM 
model is more reliable than the BPM one for all 
radial orders considered (n r = 10 to 30), being 
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Fig. 7. — Comparison of GOLF data with those of 
the TPM (crosses) and the BPM (open circles) models 
for I = 2,3,4; sizes of circles and crosses are inversely 
proportional to I. The Relative Error value is given by 

a MODEL— a GOLF x 100% 
&GOLF 



noticeable at the higher ones. 
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